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ABSTRACT 


A  method  of  calculating  the  wind  tunnel  interference  induced  by 
arbitrary  wall  configurations  has  been  developed  wherein  a  linearized 
equation  of  the  perturbation  velocity  potential  is  used  to  describe  the 
flow  in  the  tunnel.  A  series  solution  of  the  interference  velocity  poten¬ 
tial  satisfies  the  differential  equation  exactly.  The  constant  coefficients 
in  the  series  solution  are  determined  by  the  point  matching  technique 
along  the  boundary.  Results  are  presented  for  slotted  tunnel  walls  with 
homogeneous  boundary  conditions.  Solutions  presented  show  excellent 
agreement  with  the  results  obtained  by  other  methods  for  a  stationary 
as  well  as  an  oscillating  wing.  Finally,  the  method  is  applied  to  Com¬ 
pute  a  set  of  zero  interference  configurations  for  a  combination  of  slot 
openings  of  side  walls,  roof,  and  floor. 
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SECTION  I 
INTRODUCTION 


Several  investigations  have  been  made  theoretically  and  experi¬ 
mentally  to  determine  an  optimum  test  section  of  a  Vertical/ Short 
Takeoff  and  Landing  (V/STOL)  wind  tunnel  (see  references  cited  in 
Ref.  1).  Although  a  zero  interference  configuration  has  not  yet  been 
determined,  all  results  indicate  that  the  wind  tunnel  with  the  mixed 
boundaries  is  a  promising  configuration  to  alleviate  the  interference 
effects. 

As  the  first  step  to  obtain  a  theoretical  solution  to  the  interference 
problem,  a  method  is  required  to  calculate  the  tunnel  interference  on 
the  horizontal  trailing  wake  which  is  the  limiting  case  of  a  V/STOL  model. 
Unfortunately,  the  available  theoretical  methods  are  inadequate  for  the 
mixed  boundary  problem.  The  classical  image  method  has  no  proper 
image  system  satisfying  the  boundary  conditions.  In  recent  computer  in¬ 
vestigations  (Refs.  2  and  3),  a  vortex  lattice  network  has  been  utilized 
to  represent  the  model  and  wind  tunnel  walls.  However,  the  technique  is 
difficult  to  use  because  the  choice  of  the  size  of  networks  and  the  distribu¬ 
tion  of  networks  seem  ambiguous. 

In  this  report,  the  point  matching  technique  (Ref.  4)  is  introduced  to 
calculate  the  interference.  The  attractive  feature  of  the  method  is  that 
it  is  applicable  to  a  tunnel  with  an  arbitrary  cross-section  and  wall  con¬ 
figuration.  Calculations  were  made  for  the  case  of  a  horizontal  wake 
(parallel  to  the  tunnel  centerline)  and  compared  to  solutions  obtained  by 
other  methods.  In  addition,  the  wall  interference  is  calculated  at  a  point 
in  the  downstream  wake  of  an  oscillating  wing.  Solutions  are  also  ob¬ 
tained  for  the  wall  interference  for  the  horizontal  wake  case  in  a  wind 
tunnel  with  various  slot  configurations.  The  results  indicate  the  proper 
combination  of  slotted  tunnel  walls  that  will  result  in  zero  lift  interference 
for  the  case  of  a  horizontal  wing  wake. 


SECTION  II 
GENERAL  ANALYSIS 


The  flow  field,  in  terms  of  the  perturbation  velocity  potential,  is 
governed  by  Laplace's  equation  fox'  the  problem  of  the  wall  interference 
in  a  wind  tunnel.  In  the  pi'esent  analysis  a  wing  of  span  "s"  is  considered 
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and  is  represented  by  a  single  horseshoe  vortex.  The  perturbation 
potential  is  written  as: 

4>  =  4>m*  +  <£i 

where 

<£m'  =  the  potential  caused  by  the  wing  in  the  absence  of  tunnel  walls, 

4>\  —  the  interference  potential  induced  by  the  walls. 


The  potential  at  the  point  (x,  y,  z),  as  shown  in  Fig.  1,  Appendix  I, 
because  of  the  wing  at  the  origin  with  the  circulation  r  about  the  wing, 
is 


*p  m  * 


Ts  z  ^  X 

47r  y1  +  z2  (x2  +  y1  +  z2)4 


(2) 


The  potential  for  the  wake  far  downstream  from  the  wing  is 

‘Pm  =  Hltl  (f>m*  =  —  — — ^3) 

x-*m  2rr  yJ  +■  z 

which  is  independent  of  x.  This  will  be  used  as  the  basic  flow  far  down¬ 
stream  when  no  walls  are  present. 


2.1  DIFFERENTIAL  EQUATION  OF  THE  INTERFERENCE  POTENTIAL 


The  governing  equation  for  the  velocity  potential  in  incompressible 
flow  is 


d2<p  d2<b  PPtb 

<9x2  +  ~d?  + 


(4) 


If  the  tunnel  geometry  is  independent  of  x,  then  it  may  be  assumed  by 
observing  Eq.  (3)  that  the  potential  of  the  far  downstream  wake  for 
steady  flow  is 

4>  =  4>  (y»z)  (5) 


Substituting  Eq.  (5)  into  Eq.  (4),  the  governing  equation  results  in 

+  o 

<V  t?z3 

or  in  polar  coordinates,  p  =  (y2  +  an(j  q  _  tan"  (z/y), 

I  .  _L  ^  _  n 

dp2  p  dp  p2  dti2  (6) 
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The  equation  describing  the  interference  potential  is  obtained  by 
combining  Eqs.  (1),  (3),  and  (6), 

t  i_  <%+  J_  __  Q 

dp 1  p  dp  p2  dO3 


2.2  BOUNDARY  CONDITIONS 

It  is  assumed  that  along  the  tunnel  wall  there  are  longitudinal  slots. 
The  condition  to  be  satisfied  on  the  solid  portion  of  the  boundary  is  that 
there  is  no  flow  through  the  wall.  This  can  be  expressed  as 

—  =  0  at  the  solid  wall  portion  (8) 

dir 

where  n  is  the  normal  at  the  wall.  Along  the  slotted  portion  the 
pressure  is  constant  which  corresponds  to 

<$>  =  0  at  the  slotted  portion.  (9) 


Instead  of  separate  conditions  on  the  solid  and  slotted  portions  of  the 
tunnel  (discrete  slots),  an  equivalent  homogeneous  boundary  condition 
(Ref.  5)  can  be  introduced  to  satisfy  all  points  of  a  uniformly  spaced 
slotted  wall.  The  condition,  based  on  a  periodic  slot  of  width  a  and 
spacing  I,  is  as  follows: 


4>  + 


0  at  the  boundaries 


(10) 


where  K  is  the  geometric  slot  parameter 


K  -  ? '» [“*  l] 

and  a/i  is  the  open  area  ratio  as  shown  in  Fig.  2. 


(11) 


2,3  METHOD  OF  SOLUTION 

The  general  solution  of  the  interference  potential  of  Eq.  (7)  may  be 
written  in  series  form  with  undetermined  constant  coefficients,  i.  e. , 

<£i  -  jcin  sin  nO  -r  C21,  cos  ndj  {c3n/Pn  +  C4np-nj  (12) 

Since  <}>i  must  be  finite  at  p  =  0,  C^n  must  be  zero  and  Eq.  (12)  becomes 

d>i  =  J~B lnpn  sin  n0  +  B2 np"  cosndl  (13) 
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Equation  (13)  is  an  exact  solution  to  the  differential  equation  and  it 
may  be  applied  by  the  point  matching  method  to  approximate  the  boundary 
conditions  so  as  to  determine  the  constants.  If  the  series,  Eq.  (13),  is 
truncated  at  n  =  N,  then  there  are  2I\'  +  1  unknown  constants.  For  an 
arbitrary  cross-sectional  shape,  2N  +  1  points  may  be  chosen  along  the 
boundary.  Evaluation  of  either  Eqs.  (8)  and  (9)  or  Eq.  (11)  at  these  points 
results  in  2N  +  1  simultaneous  linear  algebraic  equations  for  the  undeter¬ 
mined  constants  B^n  and  B2n-  The  inversion  of  this  set  of  simultaneous 
equations  completes  the  solution  of  the  problem. 


The  mathematical  expressions  and  method  of  solution  proceeds  as 
follows:  Let  the  2N  +  1  linear  equations  be  written  in  matrix  form 


IMJ  IB  I  =  (Df 


(14) 


where 


Since 


[m]  is  the  coefficient  matrix, 

{B/  is  the  unknown  vector  and 

{Dj  is  the  right-hand  side  column  matrix 
of  this  set  of  linear  equations. 

[M]  is  a  square  matrix,  the  unknown  constants  Bn  may  be 
evaluated  as 


IB!  =  [Mr1  )D! 


(15) 


Alternately,  a  more  accurate  result  may  be  obtained  (Kef.  6)  by  the 
method  of  least  squares  and  by  choosing  more  matching  points  along  the 
boundary  than  the  number  of  undetermined  constants.  Let  the  error  on 
the  boundary  be  expressed  as 


*b  ~  (16) 

where  is  the  exact  expression  for  interference  potential  on  the 

boundary.  The  error,  eg,  is  expressed  in  terms  of  the  undetermined 
constants  Bj,  j  =  1,  2,  ...  2N  +  1,  where  (2N  +  1)  is  less  than  the  number 
of  matching  points,  Q.  The  sum  of  the  squares  of  the  errors 


S(Bp 


Q 

I 

i=  l 


eBl 


Mi, 8, 


(17) 


is  then  minimized  through  the  (2N  +  1)  relationships 

||  =  0,  k  =  1,  2,  .  .  .  2N  +  1 


(18) 
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The  resulting  system  of  algebraic  equations  is 

Q  r2N-n  -i  q 
Ji  MuBi]  =  MikOi 

[M]t  [M]  <B|  =  [M]t1D|  (19) 

The  solution  is  then  given  by  the  inversion  of  Eq.  (19). 

The  error  in  the  domain  of  the  solutions  for  the  above  two  methods 
is  bounded  by  the  boundary  error,  eg.  It  is  known  from  the  maximum- 
minimum  theory  of  harmonic  functions  that  the  absolute  value  of  error 
has  its  maximum  value  on  the  boundary.  Therefore,  the  smaller  the 
boundary  error,  the  smaller  the  error  in  the  interior  domain  becomes. 
In  other  words,  the  greater  the  number  of  points,  the  more  accurate  is 
the  solution, 


2.4  LIFT  INTERFERENCE  FACTOR 


Once  the  interference  potential  is  obtained,  the  induced  velocity 
caused  by  the  walls  can  be  calculated  at  any  point  in  the  tunnel.  The 
wall  interference  on  lift  is  expressed  as  an  angle  whose  tangent  is  the 
vertical  component  of  interference  velocity  divided  by  the  free- stream 
speed,  AW/U,  The  classical  interference  factor  8  is  defined  by 


where 


SSCl/C  =  Aaj  =  tan  Aaj  =  AW/U 

S  =  wing  area 
C]_  =  wing  lilt  coefficient 
C  =  wind  tunnel  cross-sectional  area 


The  factor  may  also  be  expressed  in  terms  of  wing  circulation  and  wing 
span  as 

8  -  CAW/2TS 


The  interference  factor  at  the  plane  of  wing  60  is 

80  =  CAW«/4rs  (20) 

where 


z=Q 


AWro  is  the  interference  velocity  far  downstream  from  the  wing  and  by 
symmetry  is  twice  the  value  at  the  plane  of  the  wing. 
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SECTION  III 
RESULTS 

3.1  RECTANGULAR  TUNNEL  WITH  SOLID  SIDE  WALLS  AND  SLOTTED 
ROOF  AND  FLOOR 

The  first  example  chosen  to  verify  the  present  method  is  the  wind 
tunnel  interference  of  a  rectangular  cross  section  with  solid  side  walls 
and  slotted  roof  and  floor.  This  case  has  been  investigated  by  analyti¬ 
cal  methods  as  well  as  an  experimental  method  utilizing  an  electrical 
analogy  to  obtain  theoretical  solutions. 

The  coordinate  system  of  the  test  section  is  shown  in  Fig.  2.  For 
the  equivalent  homogeneous  boundary  conditions,  the  overall  effect  of 
the  slotted  roof  or  floor  is  given  by  Eq.  (10), 

<6j  +  Kr  =  -  ^>m  +  Kr  at  z  =±h  (21) 

where  Kr  is  the  geometric  slot  parameter  of  the  roof  and  floor  defined  in 
Eq.  (11).  For  the  solid  side  walls,  the  normal  velocity  at  the  boundary 
is  zero. 

sT"  "  ‘ST*1’'  '  (22) 


The  wing  is  located  at  the  center  of  the  tunnel  for  convenience.  How¬ 
ever,  it  should  be  noted  that  the  method  is  applicable  for  a  wing  located 
at  any  position  in  the  tunnel.  The  potential  of  the  wing  from  Eq.  (3)  is 


=  7rT^  (23) 

Let  0  -  0  be  in  the  plane  of  the  wing.  Then,  because  of  the  symmetrical 
location  of  the  wing  and  the  asymmetrical  potential  of  the  wing  about 
0=0,  it  is  apparent  that  <f>^  (0)  =  </>i  {it  -  9 )  and  <j>^  (-0)  =  -  <^(0).  Thus  the 
solution,  Eq.  (13),  may  be  written  in  the  form 


<t>\  =  .  2 


rti  — 1#3(5 


Am  pm  sin 


(24) 


By  applying  the  boundary  conditions  of  Eqs.  (21)  and  (22)  to  Eqs.  (23)  and 
(24),  the  following  sets  of  equations  are  obtained: 


c=U3.S  Am  P° 


sin  m 0r  +  Kr  m  /s™  1  cos  (ra-1)  9r 


(25) 
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where 


and 


where 


% 

pi  =  (yj  +  h1)  8,  =  tan  1  (h/y) 


^  ,  Am  in  sin(m  -  i)  6, 

111  —  1  |  g  |  5 


r  s  2?b 

2it  p„* 


i/ 

Pxv  =  (bs  +  zJ)  =  tan"1  U/b) 


(26) 

(27) 

(28) 


Tne  constant  coefficients  Am  are  determined  by  the  point  matching 
technique  from  Eqs.  (25)  and  (27). 


The  lift  interference  factor  at  the  plane  of  wing  only  involves  the 
first  constant  Aj,  as  seen  from  Eqs,  (20)  and  (24) 


So 


CA, 

4Ts 


(29) 


The  interference  factor  is  tabulated  in  Table  I  (Appendix  II)  corresponding 
to  truncating  the  series  at  various  numbers  of  terms  to  indicate  the  con¬ 
vergence  of  the  series.  It  can  be  seen  that  there  are  only  small  differ¬ 
ences  in  the  interference  factor  as  the  number  of  terms  increases.  Table 
II  gives  the  results  of  the  series  truncated  after  the  ninth  term  for  various 
numbers  of  matching  points.  The  value  of  the  interference  factor  is  seen 
to  converge  as  the  number  of  matching  points  increases. 


The  convergence  of  the  interference  factor  for  various  height- to -width 
ratios  was  also  examined.  It  may  be  seen  from  Table  III  that  the  case 
of  X  =  1.0  (square  tunnel)  converges  most  rapidly.  It  was  pointed  out  in 
Ref.  6  that  this  may  be  partially  because  the  ratio  of  maximum  to  mini¬ 
mum  boundary  radius  (Pmax^min^  has  the  smallest  value. 

The  interference  factor  is  plotted  as  a  function  of  the  nondimensional 
slot  parameter,  Pr  =  (1  +  Kr/h)_1,  in  Fig,  3  for  various  height-to-width 
ratios.  The  value  Pr  =1.0  corresponds  to  solid  roof  and  floor,  and 
Pr  =  0.  0  corresponds  to  open  roof  and  floor.  Results  from  Refs.  5,  7, 
and  8  are  also  shown  in  the  figure.  It  can  be  seen  that  the  agreement  be¬ 
tween  the  different  methods  is  excellent. 


3.2  RECTANGULAR  TUNNEL  WITH  FOUR  WALLS  SLOTTED 

The  wind  tunnel  interference  in  the  rectangular  tunnel  with  four 
walls  slotted  is  chosen  as  the  second  example  (see  Fig.  2).  This  par¬ 
ticular  configuration  has  not  been  solved  analytically.  The  only  avail¬ 
able  solution  has  been  obtained  by  Rushton  (Ref.  8)  by  an  electrical 
analogy  technique, 
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Again,  the  equivalent  homogeneous  boundary  conditions  will  be 
used.  For  the  slotted  roof  or  floor  Eq.  (25)  is  used.  The  boundary 
condition  for  the  slotted  side  walls  is  governed  by 

di  ±  =  -  Um  ±  Kw^?\aty  =  ±b  (30) 


Combining  Eqs.  (23),  (24),  and  (30),  the  following  equation  is  obtained 


X 

m=l,  3,5 


P 

•  w 


m  3  *, 


+ 


Kw  m 


P 


m^l 

w 


sin  (m  -  1)  6 


w 


(31) 


Kw 


2  zb  - 

P  W 


where  pw  and  0W  are  polar  coordinates  as  defined  in  Eq.  (28). 


The  constant  coefficients  Am  are  determined  by  the  point  matching 
techniques  from  Eqs.  (25)  and  (31)  for  Kr  =  Kw.  The  lift  interference 
factors  are  plotted  in  Fig.  4  and  show  good  agreement  with  Rushton's 
electrical  analogy  results. 


3.3  OSCILLATING  WING  IN  A  SLOTTED  WIND  TUNNEL 


This  section  is  devoted  to  the  case  of  the  harmonically  oscillating 
wing  that  is  restricted  to  a  rotatory  oscillation  about  the  center  of  the 
chord.  The  interference  factor  in  the  infinite  wake  can  only  be  used  as 
an'  indication  of  the  amplitude,  but  not  the  phase,  of  the  interference 
factor  in  the  plane  of  the  wing.  For  a  wing  of  very  small  chord  with 
uniform  spanwise  loading,  the  expression  of  the  potential  field  as 
x— *■  ®  of  a  small  wing  oscillating  with  an  angular  frequency  u  has  been 
derived  by  Rushton  (Ref.  8)  and  Goodman  (Ref.  9)  and  is  given  by 


$ni(y»z) 


USClCO  a 
4tt  Up 


(32) 


where  Kj  is  the  modified  Bessel  function  of  the  second  kind  of  order  1. 


For  a  tunnel  of  constant  cross  section,  it  may  be  assumed  that  the 
total  perturbation  potential  of  the  far  downstream  in  a  harmonically 
oscillating  wake  is  given  by 


4>  ~- 


(33) 


Equations  (4)  and  (33)  lead  to  the  differential  equation  for  0  (y,  z) 
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or  in  polar  coordinates 


iv  +  i  +  j_  iv  =  j 

p  ap  p2  as* 


(34) 


(35) 


From  Eqs,  (1),  (32),  and  (35),  the  differential  equation  for  the  inter¬ 
ference  potential  is  given  as 


a2^,  i  a<£;  i  a1^  /&j\i 

dp2  p  dp  p2  aa2 


(36) 


The  general  solution  of  Eq.  (36)  in  series  form  becomes 


oo 

=  Mc‘  n  sin  n0  +  C 2n  cos  nSj  |C3a  In  (ap)  +  C.4nKD  (ap)j 


(37) 


where  In  and  Kn  are  the  modified  Bessel  function  of  the  first  and  second 
kind  of  order  n,  respectively,  and  a  =  ui/U. 

For  a  rectangular  wind  tunnel  with  solid  side  walls  and  slotted  roof 
and  floor,  the  boundary  conditions  are  the  same  as  that  of  the  nonoscil¬ 
lating  wing  case,  Eqs.  (21)  and  (22).  Using  the  properties  of  symmetry 
about  6  -  ir 1 2,  asymmetry  about  0-0  and  C4n  =  0,  <£■[  can  be  written  as 

oc 

0 i  -  £  A,„  Iln  [ap)  sin  ra 0  (38) 

1 1  o  f  5 

From  Eqs.  (32)  and  (38)  and  the  boundary  conditions  of  Eqs.  (21)  and 
(22),  the  constants  Am  may  be  determined  by  the  following  sets  of 
equations: 


2 

m  =  lt3,  S 


e  in  m0T  +  K 


)  c 


os  (m—1)  6r 


+  ft  [n+i  (ft)  sin  mdT  sin  9  ■ 

-  “  7T  jr  K‘  w  p*  +  Kr  [(yJ  '  hJ)  K>  W 

-  ft  h2  Ko  (ft)]j 


(39) 
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where 


and 


where 


Pr  —  (ha  +  y2)'2  et  =  tan"1  (h/y) 

/3r  -  apr 

OQ 

S  A 

m=l,  3jS 


iS. 


Im+'l  sin  m9v 


cos 


+■  m 


(j3w)  sin  (m  -  1)  0W 


USCL 

4rr 


-SjpKi  (jSw)  +  j8w  K0  (j8w 
Pw  L 


) 


=  (tf  +  z2) 


0W  =  tan  1  (z/b) 


—  aP  w 


(40) 


(41) 


(42) 


The  constant  coefficients  are  determined  from  this  system  of  equations 
by  the  point  matching  technique  along  the  boundaries  as  before. 

The  upwash  interference  at  the  centerline  in  the  infinite  wake  is 

dd>, 

AW.*  =  lim  — —  =  a  A,/2  (43) 

p  -*  0  <7z 


The  lift  interference  factor  in  the  infinite  wake  becomes 

„  -c  c  (44) 

2  USCL 

Figure  5  shows  the  interference  factor  as  a  function  of  the  reduced 
frequency  uh/U  for  various  values  of  the  nondimensional  slot  parameter, 
Pr  =  (1  +  Kr/h),  for  X  =  0.  5.  It  is  seen  that  the  interference  factor  is  of 
the  same  order  of  magnitude  as  the  nonoscillating  case  at  low  frequencies 
but  is  negligible  at  high  frequencies. 


It  is  also  interesting  to  note  the  slot  parameter  as  a  function  of  fre¬ 
quency  parameter  for  the  zero  interference  condition  as  shown  in  Fig.  6. 
It  can  be  seen  that  the  higher  the  frequency  the  greater  the  percent  open¬ 
ing  required  for  zero  interference.  For  values  of  the  slot  parameter 
less  than  0.  38,  zero  interference  is  never  reached. 
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SECTION  IV 

APPLICATION  TO  OPTIMIZATION  OF  WALL  CONFIGURATIONS 


Since  the  point  matching  technique  permits  an  interference  evalua¬ 
tion  of  wall  configurations  with  mixed  boundaries,  it  becomes  possible 
to  optimize  the  wall  configuration  of  a  tunnel  on  the  basis  of  a  number 
of  criteria.  An  example  of  applying  the  method  in  search  of  a  set  of 
zero  interference  wall  configurations  can  proceed  as  follows. 

The  first  consideration  could  be  the  tunnel  height- to- width  ratio. 

From  the  study  of  a  finite  span  wing  in  a  closed  test  section  {Ref.  3),  it 
was  found  that  the  tunnel  with  a  height-to- width  ratio,  A,  from  0.  667  to 
0.  8  gives  a  reasonably  small  rate  of  variation  of  the  correction  for  angle 
of  attack  and  velocity  along  the  lateral  and  longitudinal  axes.  In  the  in¬ 
vestigation  of  "recirculation  limits,  "  Rae  (Ref.  10)  concluded  that  a 
height-to- width  ratio  of  the  order  of  0.  8  is  probably  desirable  for  V/STOL 
testing  in  a  closed  tunnel.  Therefore  height-to- width  ratios  from  0.  5  to 
1.  0  can  be  chosen  in  a  st\idy  of  optimum  configurations. 

The  second  consideration  could  be  the  distribution  of  slot  openings  in 
the  tunnel  walls.  If  it  is  assumed  that  the  slot  parameters  of  roof  and 
floor  have  the  same  value,  Pr,  which  may  be  different  from  the  values  for 
the  walls,  Pw,  the  properties  of  symmetry  about  Q  =  irl 2  and  asymmetry 
about  6=0  for  the  interference  potential  are  preserved.  The  solution  of 
interference  potential  has  the  form  of  Eq.  (24).  The  houndary  condition 
of  the  slotted  roof  or  floor  is  given  in  Eq.  (25),  and  for  the  slotted  side 
walls  in  Eq.  (31). 

The  interference  factor  is  plotted  in  Fig.  7  as  a  function  of  the  roof 
and  floor  slot  parameter  Pr  for  A  -  0,  667  and  various  values  of  slot 
parameter  Pw.  It  can  be  seen  that  proper  combinations  of  Pr  and  Pw 
give  zero  interference.  Similarly,  the  combination  of  Pr  to  Pw  for  zero 
interference  is  obtained  for  other  height-to- width  ratios.  The  values  of 
wall  slot  parameter  that  result  in  zero  interference  are  shown  in  Fig.  8 
for  various  values  of  A .  It  is  concluded  that  the  interference  factor  is 
insensitive  to  the  opening  of  side  -walls  except  for  the  case  of  a  square 
tunnel.  This  gives  the  freedom  on  the  side  walls  to  satisfy  other  criteria. 


SECTION  V 

CONCLUDING  REMARKS 


A  new  method  has  been  proposed  for  calculating  the  wind  tunnel  inter¬ 
ference  of  arbitrary  wall  configurations  on  a  lifting  wing.  The  total  per¬ 
turbation  potential  is  split  into  two  parts;  the  potential  caused  by  the  wing 
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and  the  interference  potential  induced  by  the  walls.  The  unknown  con¬ 
stants  of  a  truncated  series  solution  for  the  interference  potential  are 
determined  by  the  point  matching  technique  along  the  mixed  boundaries. 

Calculations  of  the  wall  interference  on  steady  and  oscillating  wings 
are  made  to  evaluate  the  proposed  method.  The  agreement  is  excellent 
between  the  results  calculated  by  the  proposed  method,  the  available  data 
from  the  classical  solutions,  and  the  solutions  to  the  equations  as  ob¬ 
tained  by  electrical  analogy. 

Finally,  a  set  of  zero  interference  configurations  are  found  for  com¬ 
binations  of  slot  opening  on  the  side  walls  and  the  roof  and  floor.  The 
insensitivity  of  the  interference  to  the  opening  on  the  side  walls  leaves 
one  adjustable  variable  to  fulfill  some  other  optimizing  criterion. 


REFERENCES 


1.  Heyson,  H.  H.  "Wind  Tunnel  Wall  Effects  at  Extreme  Force 

Coefficients.  "  International  Congress  of  Subsonic 
Aeronautics,  N.  Y.  ,  A  pril  1967. 

2.  Joppa,  R.  G.  "A  Method  of  Calculating  Wind  Tunnel  Interference 

Factors  for  Tunnels  of  Arbitrary  Cross-Section.  " 

Report  67-1,  February  1967,  University  of  Washington. 

Also  NASA- CR- 845,  July  1967. 

3.  Bhateley,  I.  C.  "A  Theoretical  Optimization  of  the  Test- Section 

Geometry  for  the  Proposed  Boeing  Low- Speed  Wind  Tunnel.  " 
Report  D6- 15023  TN,  Boeing  Co.,  Airplane  Division, 

January  1967. 

4.  Conway,  H.  D.  "The  Approximate  Analysis  of  Certain  Boundary- 

Value  Problems,  "  J.  Applied  Mechanics,  27,  June  1960, 
pp.  275-277. 

5.  Davis,  D.  D  Jr.  and  Moore,  D.  "Analytical  Study  of  Blockage 

and  Lift-Interference  Corrections  for  Slotted  Tunnels  Obtained 
by  the  Substitution  of  an  Equivalent  Homogeneous  Boundary  for 
the  Discrete  Slots.  "  NACA  RM  L53E076,  June  1953. 

6  Leissa,  A.  W.  "Investigation  of  the  Utilization  of  the  Point 

Matching  Technique  for  the  Solution  of  Various  Boundary  Value 
Problems."  AFFDL-TR-66-  186. 

7.  Theodorsen,  T.  "The  Theory  of  Wind-Tunnel  Wall  Interference.  " 
NACA  Report  410,  1932. 


12 


AEDC-TR- 68-42 


8.  Rushton,  K.  R,  ’’Studies  of  Slotted  Interference  Using  an 

Electrical  Analogue.  "  ARC  R  &M  3452,  June  1965. 

9.  Goodman,  T.  R.  "The  Upwash  Correction  for  an  Oscillating 

Wing  in  a  Wing  Tunnel.  ”  J.  of  Aero.  Sciences,  Vol.  20, 
No.  6,  June  1953,  pp.  383-386. 

10.  Rae,  W.  A,  ,  Jr.  "Limits  on  Minimum  Speed  V/STOL  Wind 
Tunnel  Tests.  "  J.  of  Aircraft,  Vol.  4,  No.  3, 

'May-June  1967,  pp.  249-254. 


13 


AEDC-TR-66-42 


APPENDfXES 
J.  ILLUSTRATIONS 
II.  TABLES 


15 


h-M _  . 

H°h 

GEOMETRIC  SLOT  PARAMETER  • 


Fig.  2  Coordinate  System  and 
Longitudinal  Slots 


of  Wind  Tunnel  Cross  Section  and 


AEDC-TR-68-42 


INTERFERENCE  FACTOR 


AEDC-TR- 68-42 


19 


0 

CLOSEO 


0.2  0.4  0.6  0.6 

SLOT  PARAMETER,  Pr  AND  Pw 


1.0 

OPEN 


Fig.  4  Lift  Interference  Factor  for  Rectangular  Tunnels  with  Four  Walls  Slotted,  Pr  =  Pw 
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FREQUENCY  PARAMETER,  wh/U 

Fig.  5  Lift  Interference  Factor  os  a  Function  of  Oscillation  Frequency,  \  =  l  '2 
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Fig.  7  Lift  Interference  Factor  as  a  Function  of  the  Roof  and  Floor  Slot  Parameter,  A  =  0.667 


23 


SIDE  WALLS  SLOT  PARAMETER, 


AEDC-TR-68-42 


Fig,  8  Slot  Parometers  Required  for  Zero  Interference 
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TABLE  I 

THE  INTERFERENCE  FACTOR  FOR  VARIOUS  NUMBERS  OF  TERMS 
RETAINED  IN  THE  SERIES  SOLUTION 
SQUARE  TUNNEL  (X  =  1), 

SOLID  SIDE  WALLS,  AND  SLOTTED  ROOF  AND  FLOOR 


No.  of 

No.  of 

Slot 

Interfer ence 

Terms 

Matching  Points 

Parameter 

Factor 

P 

6 

r 

o 

9 

9 

mgm 

0.1241955 

0.0899108 

mm 

-0.1250009 

;  n 

11 

0.1240588 

0.0893592 

■8H 

-0.1249996 

■  13 

13 

mm 

0.1240902 

0.0894877 

-0.1249995 

15 

15 

0.1240832 

0.0894588 

• 

-0.1249999 

17 

17 

mgam 

0.1240847 

0.0894652 

1 

wim 

-0.1249999 

19 

19 

n 

0.1240844 

0.0894639 

11 

-0.1249999 
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TABLE  II 

THE  INTERFERENCE  FACTOR  FOR  VARIOUS  NUMBERS  OF  MATCHING  POINTS, 
SERIES  SOLUTION  RETAINED  TO  NINTH  TERMS 
SQUARE  TUNNEL  (\  =  1), 

SOLID  SIDE  WALLS,  AND  SLOTTED  ROOF  AND  FLOOR 


No.  of 

Matching  Points 
(0  £  0  £  tt/2 

Slot 

Parameter 

P 

r 

Interference 

Factor 

6 

o 

A6  * 
o 

9 

mmm 

0.1241955 

_ 

0.0899108 

— 

HS 

-0.1250009 

- 

11 

0.05 

0.1241300 

.0000655 

0.5 

0.0896548 

.0002560 

1.0 

-0.1249965 

-.0000044 

13 

■931 

0.1241119 

.0000181 

0.0895661 

.0000887 

-0.1249984 

.0000019 

15 

0.05 

0.1241063 

.0000056 

0.0895332 

.0000329 

1.0 

-0.1249994 

.0000010 

17 

0.05 

0.1241048 

.0000015 

0.5 

0.0895170 

.0000162 

1.0 

-0.1249999 

.0000005 

19 

0.1241046 

.0000002 

0.0895073 

.0000097 

■&H 

-0.1250002 

.0000003 

* 

A6  increment  of  the  interference  factor  due  to 
o 

increasing  number  of  matching  points. 
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TABLE  III 

THE  INTERFERENCE  FACTOR  FOR  VARIOUS  HEIGHT-TO-WIDTH  RATIOS 

Roof  and  Floor  Slot  Parameter.  P  =0.5 

’  v 

Solid  Side  Walls 


Size  of 
Coefficient 
Matrix* 

Height-Width  Ratio,  X 

0.5 

1.0 

9x9 

-.0398929 

0.0899108 

19  x  9 

-.0410906 

0.0895073 

□  /p  . 

max  runn 

2.236 

1.414 

All  Walls  Slotted,  P,  -  Pw  -  0.3 


Size  of 
Coefficient 
Matrix* 

Height-Width  Ratio,  X 

0.5 

1.0 

1.6 

9x9 

0.5235 

0.0430 

0.1178 

19  x  9 

0.0320 

0.0429 

0.0853 

29  x  9 

0.0322 

0.0429 

0.0853 

D  /p 

max  Kmin 

2.236 

1.414 

1.881 

The  size  of  coefficient  matrix  Q  x  N  represents 
Q  matching  points  and  N  terms  in  the  series. 
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